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For a spinor gas, i.e., a mixture of identical particles with several internal 

degrees of freedom, we derive the partition function in terms of the Feynman- 

Kac functionals of polarized components. As an example we study a spin-1 

Bose gas with the spins subjected to an external magnetic field and confined 

by a parabolic potential. From the analysis of the free energy for a finite 

number of particles, we find that the specific heat of this ideal spinor gas 

as a function of temperature has two maxima: one is related to a Schottky 

anomaly, due to the lifting of the spin degeneracy by the external field, the 

other maximum is the signature of Bose-Einstein condensation. 
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I. INTRODUCTION 

In this paper we extend our approach, introduced in for a class of interacting 
polarized quantum systems, to systems with internal degrees of freedom. After a review 
and a generalization of our methods to the unpolarized case, their feasibility is tested on 
a relatively simple confined model system of spin-1 bosons in an external magnetic field. 
The partition function, the specific heat and the susceptibility of this trapped spinor gas are 
evaluated for 100 and 1000 particles. Explicit analytical results in closed form were derived 
for 6 particles using symbolic algebra. 



In the thermodynamic hmit (A^ oo with a fixed density p = N/V), Ginibre [Q] 
provided a Feynman-Kac functional description of the quantum statistical equilibrium for 
systems with internal degrees of freedom, based on the grand canonical ensemble 0]. This 
approach, applied to quantum systems |^,^ and quantum plasma's has recently been 
generalized by Cornu to mixtures of identical particles in an external magnetic field. 
A review of the methodology can be found in |10| with emphasis on low-density Coulomb 
systems. 

The approach which we present here is based on a different description although 



both Ginibre's and our approach have in common the replacement of second quantization 
by a direct use of the permutation symmetry in a path integral. But there are essential 
differences in the way that the finite number of particles is incorporated in the theory. 
Furthermore, we treat here an unpolarized mixture of identical particles. 

Our approach avoids the thermodynamical limit and keeps the relations between the 
density of states, the partition function for particles and their generating function math- 
ematically exact ||T3[. These statistical quantities are shown to be transforms of each other, 
i.e., they are not separately defined on the basis of properties of an ensemble. The partition 
function Zj^ (/?) is the Laplace transform of the distribution of states Qn' (E) with /? adjoint 
to the energy E. The generating function S {f3, 7) of the partition function is a Z-transform 
of Zjv with the fugacity 7 adjoint to the number of particles. For a harmonic model with 
interactions 0,|13[ we could show, using the appropriate saddle-point methods reliable for the 



singular structure of the generating function at large A^, that these transforms give the well 
known relation between the fugacity -and thereby the chemical potential- and the number 
of particles. Our derivation was largely inspired by the inversion method used by Fowler 
and Darwin |[T^ to obtain the relation between the temperature and the internal energy. 
The temperature as a measure for the internal energy and the chemical potential as a char- 
acterization of the number of particles are quantities obtained from the theory, whereas in 
an ensemble-based approach (valid in the thermodynamic limit), their definitions belong to 
the theory. A similar statistical methodology has been recently proposed by Bormann et al. 



as a new approach for systems with a finite number of particles, although their method 
is only applicable for a partition function of a polarized system without any particle-particle 
interaction. 

Some models of harmonically confined systems with harmonic two-body interactions can 



be solved exactly. They have been studied using other techniques by many authors ITB-ITU 



The classical version of the model was already studied by Newton pO|,pi|. The reason for the 
exact solvability is the reduction to a center-of-mass evolution and an evolution of the other 
degrees of freedom relative to the center of mass, the latter being independent of each other 
||2^ . Despite the conceptual simplicity of this approach, the actual calculation turns out to be 
quite involved. Indeed, a general two-body interaction requires a cumulant expansion of the 
exponential in the Feynman-Kac functional with the implication that n-point correlation 
functions have to be calculated for each higher order cumulant. These point correlation 
functions have their own generating functions, which have to be inverted separately. A loop, 
as is well known from a field theoretic approach to the many-body problem, manifests itself 
in the irreducible part of a n-point correlation function. In principle, these loops can be used 
to simulate a many-body system ^ or are used in a Mayer cluster expansion [|7|,pB[ of the 
statistical quantities relevant for stability studies of low-density Coulomb systems. Although 



this prescription is general, we have applied it only to the first order cumulant ||2J], in using 
the Jensen-Feynman variational principle to optimize the free energy for realistic interaction 
potentials. The second cumulant would already require a four-point correlation function. 

Another extension of Ginibre's approach 0,0, considered in this paper, concerns the 
spin degrees of freedom. We have to project the full A^-body propagator for distinguishable 
particles on the irreducible representations of the permutation group. In the calculation 
of the partition function of the model we only need the diagonal part with respect to the 
spin degrees of freedom. This means that only those propagations have to be considered 
where particles start and end in the same spin state after some Euclidean time lapse. This 
simplification cannot be made in the calculation of a n-point correlation function. In that 
case the evolution of a spin state should be described by a continuous-time Markov process 



with discrete states |25|. These states serve as indices on the sample paths in the Feynman- 



Kac functional. Also when charges are present, the influence of the magnetic field on the 



charged particles can be taken into account with the same methods ^6 . 

The care given to the A^-body aspects in our formulation originates from the fact that 
a crossover from a density dependent behavior to a behavior dependent on the number of 
particles is experimentally accessible. It has become possible to produce assemblies of atoms 
whose internal energy is so small, that quantum effects such as Bose-Einstein condensation 
or Fermi-Dirac degeneracy become observable |p7|-0O[]. The theoretical challenge posed by 
these systems is their particular kind of confinement and the relatively low number of atoms 
involved. This low number of atoms suggests that there is a regime where the number 
of particles is more important than their density. The study of such a crossover excludes 
the thermodynamical limit as an investigation tool ||31|. The energy levels are primarily 



determined by the trapping potential, rather than by the confining volume used in most 
quantum-statistical studies. In the earliest experimental realizations of Bose-Einstein con- 
densation the system was restricted to a specific internal degree of freedom; in present day 
studies this condition has been relaxed leading to the so-called trapped spinor quantum 
gases P^^^, where an equilibrium over the internal degrees of freedom can be reached. 

In the present paper we illustrate the calculation techniques for a boson gas with three 
internal degrees of freedom. A preliminary report on the results of this investigation can 



be found in The paper is organized as follows. In section II we derive the quantum 
mechanical partition function of a mixture of identical particles, given only the total number 
of particles. In order to achieve this goal, we go through a number of steps: first we derive a 
partition function conditioned on a particular state of the internal degrees of freedom and a 



given number of particles as a Feynman-Kac functional p6|-p8| that can be obtained explic- 
itly by path integral calculations, at least for the model under consideration. The resulting 
propagator is projected on the symmetric or antisymmetric irreducible representation of the 
permutation group according to the boson or the fermion character of the particles. The 
cyclic decomposition of the permutations imposes a constraint on the summations over the 



cycle lengths. This difficulty is circumvented by introducing a generating function for each 
conditioned partition function, which can be inverted. The partition function of the mixture 
is then obtained as an appropriate combination of the conditioned partition functions, from 
which the free energy and the other thermodynamical quantities of the mixture result. We 
describe these steps using an harmonic boson model, supplied with an homogeneous external 
magnetic field, i.e., with a Zeeman splitting of the spin degrees of freedom. In section III the 
free energy, the susceptibility and the specific heat will be obtained with emphasis on the 
low temperature properties and the thermal fluctuations of a boson gas. Finally, in section 
IV we discuss the method and the obtained results, and we conclude the paper with a brief 
summary of the work. 



In this section we provide the necessary background material to formulate a quantum 
statistical theory of mixtures of a fixed number of identical particles with internal degrees 
of freedom. Using the same approach as we developed |jl|J^ for particles without internal 
degrees of freedom, we have to describe the state space and the evolution on that state space 
using a process and a Feynman-Kac functional averaged over that process. We also have to 
indicate how the projection on the irreducible representations of the permutation group has 
to be done in order to ensure the indistinguishability of those particles that are in the same 
state of their internal degrees of freedom. Finally, we give an expression for the partition 
function of the mixture. 



The quantum model which we consider contains identical particles kept together in a 
confining potential given by 



II. QUANTUM STATISTICS OF MIXTURES 



A. The propagator of the harmonically confined model 




(1) 



m=—s 



k=l 
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where rm,k is the position vector of the k^^ particle in state m. (Natural units with h and 
the particle mass equal to unity are used througout this paper.) The frequencies are 
related to the curvature of the parabolic confinement potential, which in principle might be 
different for each internal degree of freedom. We assume that there are N^n particles that 
occupy the state m. The total number of particles is obtained from 



s 



N^J2^rn. (2) 

m=—s 

In the absence of interparticle interactions and for distinguishable particles, the propagator 
for the spatial degrees of freedom can be written as a product 

Kn{{r':^,k},f3\{r'^,,},0)= JJ JJ K {r^^,, (3K^„0)\^^ , (3) 

m=—s fc=l 

where {r[^ ^^j represents a configuration of all the particles in the internal state m, and with 

w {rl + rl) cosh - 2ra ■ n 
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w \ ^ 

K{ra,f3\n,Q)\.= ( ^_^^^^^^_ ] exp 



sinh w(3 



(4) 



2n smhwP 

We assume that the internal degrees of freedom behave like a spin in an external homo- 
geneous magnetic field B, described by the Hamiltonian 

n = -/xB • S. (5) 

Introducing for each particle the states characterized by the value m of the spin component 
along the quantization axis 

SzlsjUi) — mlsjUi) , m — —s, —s + 1, . . . , s, (6) 

one finds that the Euclidean time evolution for the spin states of the particle is given by: 

ICs (m", /3\m', 0) = (s, m"\ exp -pH \s, m') . (7) 

This propagator can be written in the spectral representation by introducing a unitary 
matrix U that diagonalizes H, 

unu-^ ^-l^\B\S, , (8) 



as follows 

s 

Ks (m", ^|m', 0) = ^ (s, m"\ [/"^ |s, m) e^""^ {s, m\U\s, m') , (9) 

m=—s 

with uJm = fJ'fn |B| . 

Disregarding the statistics for the time being, the propagator for all the (distinguishable) 
particles including the internal degrees of freedom thus becomes: 

r K (v" . ii\-r' . nV 



^D({r':,J,/3|{4,J,0)= n 



(10) 



k=l 

The quantum evolution of the system is now defined as follows: the particles are confined 
by a harmonic potential that may differ according to their spin state m, and the state m of 
the j*'' particle evolves as a spin-m state in a magnetic field. This model is an idealization 
because interaction between the spins of different particles is left out, spin-orbit coupling 
or pair formation is not considered, and there is no spatial two-body interaction. But in 
the assumption that these interactions can be included, this simplified model still has to be 
studied as a zero order approximation. 

B. The permutation symmetry and the generating function 

From the propagator Kd for the distinguishable particles, the propagator K^ixt of the 
mixture, taking into account the fermion or boson statistics, can be obtained by using the 
antisymmetric or symmetric projection as documented in [0,|l^: 



Krmxt i{r'U},/3\ {4, J ,0) = n e^-^--K, ({C^J ,/?| ,0) . (11) 

m=—s 

The expression for the propagator Kj (|r^ j^.} , P\ {r'^ k} 5 O) of identical particles in the same 
spin state m is 

{K,k} , P\ K,k} ' 0) = ^ E n ^ ( VP(fc)' Pl^kk, 0) In ' (12) 

^^'"^ p k=i 
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which represents a Feynman-Kac functional of the polarized components with spin state m, 
including all their permutations P. Knowing the propagator, the partition function for a 
given distribution N^s, ■ ■ ■ ,Ns can be expressed as a multiple integral: 



Z{P\N.„...,N,) 



s N„ 



n n * 



r„,» TT e''""'""'K;({r„,i.},/J|{r„,i.},0). (13) 

_m=—s k=l ' 

The calculation for each spin degree of freedom then proceeds in complete analogy with the 
polarized case leading to a partition function of a mixture with a given composition of 
internal degrees of freedom 



(14) 



where Z (/?|A^m) is the partition function of spin-polarized particles: 



Z(d\N„ 



'Nrr 



.k=l 



n 



({r„,fc} {rm,fc} ,0) 



(15) 



Of course, if the total number of particles is fixed, all combinations of the A^_s, . . . ,Ns 
(with J2m=-s ^rn = N) are possible and we obtain the following expression for the partition 
function of the mixture: 

N N N 



3rm.t{m)= Yl Yl ••• 5^ C'(iV_„...,iV,)Z(/5|iV_„...,iV,), 

N_s=ON_s+l=0 Ns=0 



C(iV_„...,iV,) 



(16) 



(17) 



N_JN_,+,\...NJ ' 

where C {N_s, . . . , A^^) is a multinomial coefficient in which B (A^ = Ylm=-s expresses 
the constraint on the summations, i.e., (x) = 1 if the logical variable x is true, and 
(x) = if the logical variable x is false. Using the expression for the joint partition function 
"bmixt in terms of the polarized partition functions Z {P\Nm) one readily obtains: 



N N 



N 



imixt 



:i8) 



C{N_s,...,Ns) n Z(/5|Arje^™^-" 

This is the partition function that will be analyzed in the next section. We consider 
there a gas of atoms with Bose-Einstein statistics confined in a harmonic potential, without 
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interparticle interactions, and with the abihty to adjust their internal degrees of freedom in 
order to attain thermodynamical equihbrium. If the system is prepared in such a way that 
the total z component of all the spins can be fixed, then this constraint can be taken into 



account using a technique developed in |^. In this case one has also to calculate expression 
(|I|) first. 

III. A SPIN-1 BOSE GAS 

In this section we illustrate the feasibility of the proposed approach on a set of spin-1 
bosons confined in a harmonic potential, based on the partition function (0) of the mixture. 
For simplicity, the many-body interaction, studied in to allow for oscillations of the 
center of mass with a different frequency than that of the oscillations of the internal degrees 
of freedom, is not taken into in consideration here. This leads to a simpler albeit non trivial 
model. The calculation techniques can almost entirely be based on those that have been 
used for the polarized case 



The generating function S (/?,«) = X]Xi=o ^ (/^l^*") "^^^^ ^'^^ partition function 
Z(/?|A^m) of each component in the mixture is known from these previous calculations, 
and given by: 

1^7(^3^) with& = e-^^-. (19) 

It should be noted that u and h depend on the parameters of the spin-polarized subsystem. 
In the actual calculation we consider the same confining frequencies = w for each 
subsystem (not exploiting the possibility of different confining potentials for different spin 
states). 

There exist many ways to derive recursion relations between the partition functions with 
a different number of particles A binomial expansion of the generating function 



almost immediately leads to 

^' ' ' N 

n=0 



1 / 5(^-")/2 \ 3 



Combining the solutions for Z(/3|A^) from the recursion relations, the partition function of 
the mixture becomes: 

(21) 

which is the key result of the present paper. Once this function is known the thermodynam- 
ical quantities can be obtained in a straightforward manner. 



A. Specific heat and magnetic susceptibility 

The numerical representation of the results has to cope with very large numbers. For 

accuracy reasons the following transformation proves to be useful: 

'hrm.t (/5|iV) = iV!e^'^-^6i^5^,,„ (22) 

with 

TV iV-m 

Smixt = exp {-^Us {m + In)) (23) 

^-^ ^-^ IN — m — ny.minl 

where the scaling of the polarized partition function is given by: 

Z(/5|iV) = 6i^Xiv, (24) 
with the appropriate recursion relation: 

YAT = — >^ — o, with Yn = 1. (25) 

m=0 V 

Remembering the definition ujm = fJ^m |B| , it is clear that ujg denotes the maximal possible 
frequency due to the magnetic field. 

The expected number of particles {Nm {P\N)) in each spin state m then takes the fol- 
lowing form suitable for numerical computation: 
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^ TV N-m 

"^'^^^ m=0 n=0 ^ ' 

^ N N-m 

i'^MN)) = ^E E -nexp(-/5.4n» + 2n)) (27) 

m=0 n=0 ^ ' 

(iV_i (/3|iV)) =N-{No im)) - {Ni im)) ■ (28) 
The same parameterization can be used to obtain the internal energy of the mixture 

Urm:ct{m,tJs) ^ -^ln{3rm.t{m)) , (29) 

leading to 

V2 y Smixt^^^^Q [N - m - ny. m\ n\ 

—m—n 

+ Urn + Un + uJs{m + 2n)) (30) 

(with Um the internal energy of the particles in spin state m) . The specific heat 

C{P\N,u;s) = ■^Umi,t{m,^s) (31) 

and the magnetic susceptibility -^U^nixt : proportional to '^^^^'^"■^ , are easily obtained 

from this expression. 

In fig. 1 we show the specific heat, the magnetic excess specific heat (i.e. the contribution 
to the specific heat due to the magnetic field) and the magnetic susceptibility for 100 bosons. 
For 1000 bosons the same quantities are shown in fig 2. The temperature is expressed in 

/ \ 1/3 

units of Tc = ^ ( ' ^^^^^ ^ frequency parameter of the confinement, k is the 

Boltzmann constant and C (3) = 1.202 056 903 is a Riemann zeta- function. The frequency 
parameter for the internal degrees of freedom is expressed in units of w: 



Using these units the expression for b becomes b — exp ^— j (^^^ ^ with t — ^ 
exp {—f3uJs) = Because there is a substantial dependence on the magnetic field strength, 
the infiuence of the magnetic field can be illustrated by the redistribution of the particles over 

11 



the internal degrees of freedom, giving rise to the magnetic excess specific heat. The magnetic 
susceptibihty clearly illustrates the dependence of the internal energy on the magnetic field 
strength. 



B. The Schottky anomaly 

In the specific heat plotted in fig. 1 and fig. 2 we have identified the low temperature 
maxima as Schottky anomalies. This contribution to the specific heat is attributed to the 
hfting of the degeneracy of the levels of the spin degrees of freedom, due to the magnetic 
field. The effect occurs irrespective of the boson statistics, as will be illustrated below. 

For up to 6 particles we calculated all polarized partition functions by symbolic algebra. 
This allows to study the free energy and the derived quantities exactly. Because the same 
calculation can be performed using Maxwell- Boltzmann statistics, we may isolate the effect 
of the Bose-Einstein statistics on the Schottky anomaly. 

For numerically exact calculations it is useful to factorize out the dominant behavior 

u— 

Z{p\N)^z^{P) ' (32) 

[ln=l 

The recursion relation for the resulting polynomials zn {P) becomes then: 

^.OT=^E^~W 7!',L)3 ■ (33) 

The polynomials zjv {P) ior N — 1, . . . , 6 were obtained by integer arithmetic and the 
results of this calculation are listed in table 1. For low temperatures {kT -C w/2) the specific 
heat is plotted in fig. 3 as a function of the temperature and of the magnetic field. The 
contribution associated with the lifting of the degeneracy and the evolution to a single 
polarized state with increasing magnetic field clearly manifests itself. We have also studied 
the specific heat of 6 distinguishable particles with spin degrees of freedom. The result is 
shown in fig. 4, and it also exhibits a maximum due to the Schottky anomaly. Comparing 
the specific heat of distinguishable particles with that of bosons, it should be noted that the 
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Bose-Einstein statistics only weakly influences the magnitude of Schottky anomaly. But its 
relative contribution is more pronounced for Bose-Einstein statistics than it is for Maxwell- 
Boltzmann statistics because the boson character suppresses the energy fluctuations at low 
temperature leading to a smaller specific heat. 



IV. DISCUSSION AND CONCLUSION 

First we review our method, then we comment on our results for the spin-1 model and 
finally we will conclude with a brief indication of the potential use of the method. 

A first striking difference between the method worked out here and more standard ap- 
proaches to the many-body problem is that we incorporate the statistics using directly the 
representation theory of the symmetric group instead of the more common second quantiza- 
tion. As we indicated before, Ginibre's approach is based on the grand canonical ensemble, 
whereas we based our method for particle-conserving systems on the distribution of states, 
which we transformed to the partition function and its generating function. The tempera- 
ture and the chemical potential obtain their standard meaning from the saddle-point method 
used to invert the transform. The numerical analysis for the spin-polarized case published 



earlier P,p!3[| indicates that when there are sufficiently many particles, corrections to the sad- 
dle point inversion become negligible. We believe that this remains the case for unpolarized 
systems but we did not prove this yet. 

When internal degrees of freedom are involved, we argued that a Feynman-Kac func- 
tional averaged over an indexed Brownian motion describes the evolution of the system. 
The projection on the appropriate irreducible representation of the permutation group leads 
to a multinomial combination of the spin states in the partition function. The index process 
to distinguish between different spin states is not used explicitly, because only the diago- 
nal part of the spectral representation of the propagator of that process is needed for the 
thermodynamics. The fact that we can construct a process for the evolution of the internal 
degrees of freedom is nevertheless important, not only for a better understanding of the 
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behavior of the system at long Euchdean times (low temperature) but also for consider- 
ing correlations between components in different internal states. We indicated briefly that 
two-body interactions can be taken into account by a cumulant expansion of the exponen- 
tial in the Feynman-Kac functional. This implies that one-point, two-point and n-point 
correlation functions have to be calculated. Other expansions, like e.g. the Mayer cluster 
expansion, can presumably also be carried out in our formalism, but this point deserves 
further investigation. 

In order to demonstrate the feasibility of the approach, we gave an example based on an 
exactly soluble harmonically confined spin-1 boson gas in a magnetic field and calculated the 
internal energy, the magnetic susceptibility and the specific heat for the system in equilib- 
rium. Furthermore, although we illustrated the technique without two-body interactions, it 
should be noted that the variational method that we have applied to the spin-polarized case 
can also be used to study mixtures of identical particles with different degrees of freedom. 

The example itself, a noninteracting boson gas, exhibits Bose- Einstein condensation char- 
acterized by one of the maxima in the specific heat as a function of temperature. A second 
maximum is due to the Schottky contribution to the specific heat. It should be noted that 
the Schottky anomaly indicates that the z-component of the total spin, which is zero in 
the ground state, acquires values different from zero. The fact that this occurs at a lower 
temperature than the BEG for an ideal system may of experimental importance. 

We conclude that the methods which we developed for polarized systems of identical 
particles can also be extended to non-polarized systems, as shown in the present paper. 
The basic tools are a Feynman-Kac functional with an adapted process, the projection on 
the symmetric or anti-symmetric representations, and the path integral evaluation of the 
averages. 
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Table 



N 


Polynomials zn (/?) in 6 = e [see Eq. (^3|)] 



1 
2 


1 

1 

1 + 362 


3 


1 + 36^ + 7b^ + 66'' + 66^ + 106^ + 36^ 


4 


1 + 36^ + 76^ + 186^ + 216^ + 476^ + 576^ + 876^ + 806^ + 876^° + 636^^ 
+626^2 + 276^3 + 156^^ + 6^^ 


5 


1 + 36^ + 76^ + 186'' + 396^ + 746^ + 1386'^ + 2526^ + 3966^ + 6006^° + 8286^1 
+10876^2 _^ 1314513 _^ 1503614 j_ 1585515 _^ isgofe^^ + 14166^'^ + 11976^^ + 9216^^ 
+6696^0 + 4136^1 + 2316^^ + 1056^3 + 37524 _^ gj25 


6 


1 + 36^ + 76^ + 186"^ + 396^ + 996^ + 1806^ + 3906^ + 7156^ + 13236^° + 22146" 
+37136^2 _^ 56646^3 + 85896^^ + 12 1296^^ + 16 6836^'^ + 216066^'^ + 273126^*^ + 32 4426^^ 
+37 5696^0 + 41 1476^1 + 43 6746^2 + 43 8576^3 + 42 7566^4 + 39 2196^5 + 34 8546^6 
+29 1206^7 + 234366^8 + 17 565629 + 12 7336^0 + 83886^1 + 53046^2 + 29946^3 _^ 1590534 
+6876^5 + 2906^6 + 7563'^ + 156^^ 



TABLE I. In this table we list the polynomials zm {(3) for N = 1, ... ,6. The coefficients are 
obtained using integer arithmetic and symbolic algebra. 
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Figure captions 

Fig. 1: For 100 bosons the specific heat per particle (a), the magnetic excess specific heat 
per particle (b), and the magnetic susceptibility (c) are shown as a function of the 
temperature and the magnetic field. 

Fig. 2: For 1000 bosons the specific heat per particle (a), the magnetic excess specific heat 
per particle (b) , and the magnetic susceptibility (c) are shown on the same scale as in 
figl. 

Fig. 3: The specific heat per particle is shown for 6 bosons distributed over 3 spin states in 
thermal equilibrium. The maximum in the (C, yy',,)-plane with fixed T is identified as 
the Schottky anomaly. This is an entropic effect due to the lifting of the degeneracy 
of the spin degrees of freedom by an external field. 

Fig. 4: The specific heat per particle for 6 distinguishable particles is shown with Maxwell- 
Boltzmann statistics. Note the presence of the Schottky anomaly. This figure should 
be compared with figure 3 where the same quantity is shown taking the Bose Einstein 
statistics into account. 
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Figure 1 (a) 




Figure 1 (b) 




Figure 1 (c) 




Figure 2 (a) 




Figure 2 (b) 




Figure 2 (c) 




Figure 3 




Figure 4 



